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Abstract. High-quality if-band spectra of strongly reddened point sources, deeply embedded in (ultra-) compact 
H II regions, have revealed a population of 20 young massive stars showing no photospheric absorption lines, but 
sometimes strong Br7 emission. The Br7 equivalent widths occupy a wide range (from about 1 to over 100 A); 
the line widths of 100-200 km s _1 indicate a circumstellar rather than a nebular origin. The isT-band spectra 
exhibit one or more features commonly associated with massive young stellar objects (YSOs) surrounded by 
circumstellar material: a very red colour (J — K) <; 2, CO bandhead emission, hydrogen emission lines (sometimes 
doubly peaked), and Fell and/or Mgn emission lines. The large number of objects in our sample allows a more 
detailed definition and thorough investigation of the properties of the massive YSOs. In the (K, J — K) colour- 
magnitude diagram (CMD) the massive YSO candidates are located in a region delimited by the OB zero-age 
main sequence, Be stars, Herbig Ae and Be stars, and B[e] supergiants. The massive YSO distribution in the CMD 
suggests that the majority of the objects are of similar spectral type as the Herbig Be stars, but some of them are 
young O stars. The spectral properties of the observed objects do not correlate with the location in the CMD. The 
CO emission must come from a relatively dense (~ 10 10 cm~ 3 ) and hot (T ~ 2000 — 5000 K) region, sufficiently 
shielded from the intense UV radiation field of the young massive star. The hydrogen emission is produced in 
an ionised medium exposed to UV radiation. The best geometrical solution is a dense and neutral circumstellar 
disk causing the CO bandhead emission, and an ionised upper layer where the hydrogen lines are produced. We 
present arguments that the circumstellar disk is more likely a remnant of the accretion process than the result of 
rapid rotation and mass loss such as in Be/B[e] stars. 
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1. Introduction stars l|Bonnell et al.l Il998l) . Alternatively, non-spherical 

accretion through a disk could solve the "rad iation pres- 

The observational study of massive stars at the earliest gur{J problem » ( e .g. lYorke fc SonnhalterlEool . Therefore, 

evolutionary phases is a rapidly expanding field. Although thfi detection of circumst ell a r disks around young massive 

over the past decades significant progress has been made in starg ig regarded as an esS ential step in understanding the 

unravelling the formation process of low-mass stars (e.g. formation of the most massive stars _ observations at cen- 

i Shuet alj | 1987|) , it is not well understood how massive timeter and millimeter wav elengths suggest that rotating 

stars form. The contraction timescale is short, so that al- circumste llar disks a re present ar ound h igh-mass young 



ready very early in their formation massive stars will pro- stella r obiects (e g | Minier et al | IShepherd et al 

duce a copious radiation field that may hamper or even Eolnl | B eltran etll]l2004t IChini et al.ll2004t Uiang et al 



reverse the accretion process (jWolfire fc Cassinellilll987j) . ^qqeJ IPatel et alJl2005(l 

This has led to the suggestion that stars more massive 

than ~ 10M Q cannot form through (spherical) accretion Newl y born massive stars are dee ^ embedded inside 

alone, but instead form by collisions of intermediate-mass their natal cloud ' obscured from view h Y 10 - 100 ma S nl " 

tudes of visual extinction. One of the first observable sig- 



Send offprint requests to: A. Bik (abik@eso.org) natures is the radio and infrared emission produced by the 

* Based on observations collected at the European Southern ultra- (or hyper-)compact H II region that emerges when 

Observatory at Paranal, Chile (ESO programmes 64.H-0425, the contracting star becomes hot enough to ionise the sur- 

65.H-0602, 68.C-0652 and 69.C-0448). rounding material. Although radio and infrared observa- 
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tions provide information on the amount of ionising radia- 
tion, these diagnostics do not allow an accurate determina- 
tion of the photospheric properties of th e embedded young 
massive star(s) (for reviews, see e .g. Churchwell 1991; 
lHansonlll99atlGarav fc Lizandll999UChurchwel]ll2002l) . 

Our strategy aims at the direct observation of the ion- 
ising star(s) of (ultra-) compact Hli (UCHII) regions at 
near-infrared wavelengths, where the extinction is strongly 
reduced (Ak — O.lAy), but the (reprocessed) emission by 
dust at longer wavelengths has not yet set in. Recent de- 
velopments in near-infrared instrumentation have made 
this approach feasible. By first identifying the candidate 
ionising stars in IRAS sources with UCHII colours (cf. 
iKaper et al1l2006tl a nd su bsequent if-band spectroscopic 
follow-up ijBik et aJI EoOSI. their properties can be deter- 
mined. 

It turns out that our sensitivity allowed us to detect 
a bout half of the UCHIIs in the near- infra red. As shown 
inlBik et all (|2005jklKaner et alJ l(200fit) and lHanson et all 
(2002), the UCHII regions are always located in larger, 
extended HII regions. Apart from the near-infrared coun- 
terparts of some of the UCHII regions, we have spectro- 
scopically identified two other types of objects in those 
extended regions. The first type includes the photospheri- 
cally detected OB stars. Their photospheric properties are 
determined by applying the spectra l classification criteria 
as proposed bv lHanson et alJ l)l996(l and by measuring the 
shape of the line profiles. These objects cannot be distin- 
guished from normal "naked" OB stars. This suggests that 
these stars have already reached the main- sequence phase 
ijHanson et al .1120021 iBik et a.l]l2003l 120051) . 

In this paper we focus on the spectroscopic proper- 
ties of the second type of objects; 20 objects that we en- 
countered in our survey of southern UCHIIs. Their spec- 
troscopic and photometric properties differ from normal 
main-sequence OB stars and exhibit features commonly 
associated with massive young stellar objects. In our defi- 
nition, a massive YSO has already contracted to stellar di- 
mensions and likely commenced with core hydrogen burn- 
ing, thereby separating a massive YSO from a high-mass 
protostar which would still be in its contraction phase. 

The observational characteristics of massive YSOs are 
not well defined, but include: (i) a red continuum, likely 
due to hot dust, (ii) a location inside a star forming re- 
gion and (iii) an emission-line spectrum (Hi, Fell, Mgn, 
Nai, CO first overtone emission). The detection of emis- 
sion in the first overtone bands of CO (at 2.3 /xm) is an 
important signature, as it is considered to be an indication 
for the presence of a dense circumsteiiar disk. CO first- 
overtone emissi on has rarely been de tected in any astro- 
nomical object i McCregor et al.ll988h . The first detection 
of CO bandhead emission in a massive YSO was reported 
bv lScoville et all lll979l Il983h in the B e cklin-Neugeb auer 
(BN) object. iGeballe fc Perssonl (|l987|) . ICarrl (|l989|) and 
IChandler et alJ l)l993|) detected CO emission in another 
four high-luminosity YSOs : S106, NGC 2024-IR S2, S140 
and GL2789. Modelling bv lChandler et alJ l|l995h showed 
that vibrational temperatures of ~ 3500 K are required to 



produce the relative band strengths and that the size of 
the emitting region has to be quite small (~ 1 AU) and 
dense (n ~ 10 10 cm" 3 ). A circumsteiiar disk could pro- 
vide this environment, shielding the CO from the hot stel- 
lar radiation field. Alternatively, the CO might be shock- 
heated due to the interact ion of a stellar wind with the 
surrounding cloud (see also lScoville et alJll983tlTanll2004L 
in the case of the BN object), or emerging from the stellar 
wind itself l|Chandler et al.llQ95^ . It is not clear why some 
YS Os show CO emissio n and others do not. 

lHanson et alJ l|l997j) analysed the young (< 1 Myr) 
stellar population of M17 and uncovered a handful of mas- 
sive YSOs. if -band spectroscopy revealed that the hottest 
O stars are naked, while some of the later O and B stars 
show clear disk signatures (near-infrared excess, CO emis- 
sion, double- lined Pa<5 profiles). Similar conclusions were 
reached by Blum and collaborators based on their system- 
atic study o f the stellar content of obscured galactic giant 
H II r egions llBlum et al.ll999ll2000H200lHFigueredo et alJ 



11 11 IVj^lWllU \ 111 Ulll »j U Lll ■ _1_ i_/ i_l i_l . U\JVJVJ} iJWW-1.) -L J.C IUjI L 

l2002HConti Blumll2002HFigueredo et a,l.ll2005l) 

In the following we present if -band spectra of 20 can- 
didate massive YSOs deeply embedded in UCHIIs, ob- 
tained with ESO's Very Large Telescope and ISAAC. We 
study the spectral characteristics of these potentially very 
young massive stars, and search for signatures that may 
reveal information on the process of their formation. In 
the next section (Sect. |2J we introduce our sample and 
describe the observations and data reduction process. In 
Sect. [3] the if-band spectra and near-infrared photome- 
try are presented. In Sect. 01 we discuss some well-studied 
objects showing similar signatures as encountered in our 
massive YSO sample. In Sect.0we put the observed prop- 
erties of our massive YSO sample in perspective and pro- 
vide arguments that we have detected very young massive 
stars of which at least some are surrounded by a circum- 
steiiar disk, most likely remnant of the accretion process. 
We end this paper with a summary of our conclusions 
(Sect. EJ. 

2. Sample definition and data reduction process 

The candidate massive YSOs discussed in this paper 
are selected from a survey of near-infrared bright point 
source s deeply embedde d in high-mass star forming re- 
gions ijKaper et alJl200't?) . Based on their position in the 
colour-magnitude diagram, the most luminous and reddest 
sources, and thus the potentially newly born massive stars 
are identified. A medium-resolution (R — 10, 000) if-band 
spectrum has been obtained with ISAAC on the ESO VLT 
at Paranal, Chile. This spectral resolution, as well as the 
obtained high signal-to-noise ratio, permits an accurate 
correction for telluric absorption lines and enables us to 
discriminate between stellar and nebular emission. Given 
the demand on spectral resolution, a narrow slit (0.3") is 
used, thus limiting the if -band magnitude of the targets 
to brighter than K ~ 13. 

A number of objects turn out to be late-type fore- 
ground stars, but the majority are of early (OB) spectral 
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type ijBik et alj E0051. Many show a photospheric spec- 
trum for which the if-band spectral type has been deter- 
mined, resulting in a reliable assessment of their stellar 
parameters. Apparently, these stars, though embedded in 
a young Hn region, have already become "normal" OB 
main sequence stars. Several objects, however, do not show 
a photospheric spectrum, but exhibit a broad Br7 emis- 
sion line; these objects define the sample of the present 
paper (Table QJ. The photometric and spectral properties 
of the 20 candidate massive YSOs are outlined in the next 
section. 

Long-slit (120 ") if-band spectra of these objects are 
obtained in two wavelength settings, the first setting with 
a central wavelength of 2.134 /im including Br7 ("Be- 
setting"). The spectrum covered in this setting was used 
to identify the objects. To further investigate their po- 
tential massive YSO nature, for 15 objects we also ob- 
tained spectra in another setting covering the first four CO 
first-overtone bands (central wavelength 2.234 /im "CO- 
setting" ) . The presence of CO first-overtone emission is an 
important signature of massive YSOs, and modelling these 
bands provides information on the geometry and n ature 
of their circumstellar material fe.g. lBik fc Thill2004h . 

The observations in the Br7-setting were carried out 
partly in visitor mode (March 18-20, 2000) and partly in 
service mode (between May 8 and July 16, 2000). The 
observations in the CO setting were also obtained in two 
runs, one visitor mode run (February 21-23, 2002) and a 
service mode run in the period July 7 to August 2, 2002. 

The spectra are electric ghost, dark- and flat-field cor- 
rected. The wavelength calibration of the spectra taken 
in the Br7 setting is performed using arc spectra with 
the IRAF 1 task identify. For the spectra obtained in the 
CO setting the wavelengt h calibration is perform ed using 
telluric OH emission lines jRousselot et, al.ll2000l) . The ac- 
curacy of the wavelength calibration is ~ 3 km s _1 . In 
order to correct for telluric emission the object was "nod- 
ded" between two positions on the slit (A and B) such 
that the background of position A is recorded while the 
source is at position B and is subtracted from the obser- 
vation where the source is at position A, and vice-versa. 
The offset between the A and B positions was chosen such 
that the stars 2 on the slit do not overlap. In some objects 
the nebular emission (in, e.g., the Br7 line) is extended 
so that the emission in the A image is overlapping with 
the emission in the B image. After subtraction this results 
in artifacts at the position of the nebular lines. We have 
inspected the 2D spectra of all sources for the presence 
of nebular emission and conclude that in most cases the 
contribution (positive or negative) of nebular emission to 



1 IRAF is distributed by the National Optical Astronomy 
Observatories which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 

2 To increase the efficiency of our program the orientation of 
the slit was chosen such that at least two stars are observed 
simultaneously. 



the stellar (Br7) emission is negligible. However, in a few 
cases the nebular contribution is significant (~ 50 % of 
the total flux in the case of 06412nrl21, 17136nr649 and 
18006nr766, see Fig.EJ. 

Telluric absorption lines were removed using stars of 
spectral type A, observed under identic al sky condition s 
(for the list of standard stars, we refer to lBik et alJl2005|) . 
First, the only photospheric line (B1-7) in the spectrum of 
the A star needs to be divided out. It turns out that the 
best result is achieved when first the telluric features are 
removed from the if-band spectrum of the telluric stan- 
dard using a high-resolution telluric spectrum (obtained 
at NSO/Kitt Peak). This spectrum is taken under very 
different sky conditions, so a lot of remnants are still vis- 
ible in the corrected standard star spectrum, but without 
this "first-order" telluric-line correction, a proper fit of 
Br7 cannot be obtained. The Br7 line is fitted by a com- 
bination of two exponential functions. The error on the 
resulting Br7 equivalent width (EW) of our target star 
is about 5 %. After the removal of the Br7 line of the A 
star, the telluric lines are removed by taking the ratio of 
the target spectrum with the obtained telluric spectrum. 
In the CO setting no stellar features are present in the tel- 
luric standard star spectrum. The telluric-line correction 
is done using the IRAF task telluric, which allows for a 
shift in wavelength and a scaling in line strength yielding 
a more accurate fit. The task uses a cross-correlation pro- 
cedure to determine the optimal shift in wavelength and 
the scaling factor in line strength, which can be changed 
interactively. The shifts are usually a few tenths of a pixel 
(1 pixel corresponds to 16 km s _1 ); also the scaling factors 
are modest (~ 10 %). 



3. .fT-band photometric and spectral properties of 
massive YSO candidates 

In this section we present the if -band spectra of 20 ob- 
jects selected from our sample of obscured point sources in 
high-mass star-forming regions that display a broad, spec- 
troscopically resolved Br7 emission line. In Table[T]the ob- 
jects and some of their observed properties are listed. In 
the remaining part of the paper we will name the objects 
after the first 5 digits (i.e. the right ascension) of the IRAS 
point source they are associated with, together with a 
number based on our photometry (e.g. object 227 in IR AS 
06058 +2138 we refer to as 06058nr227, cf. iKaper et all 
2006). These objects are further characterised by a con- 
tinuum slope that is on average much redder than that 
of the OB stars, that are also present in most of these 
regions, indicating that these objects possess an infrared 
excess. We obtained photometric data in two near-infrared 
bands (J and K). That is not sufficient to discriminate an 
IR excess from extinction in a colour-colour diagram. 

In the following we first discuss the photometric prop- 
erties of these objects. After that, we evaluate the spec- 
tral features detected in the Br7 setting; subsequently, the 
spectra obtained in the CO setting are discussed. 
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Table 1. Sample of candidate massive YSOs. Column 1: IRAS point source catalo gue number; column 2: Name 
corresponding Hll region; column 3: Number jf -band point source l|Kaper et al . 2006); columns 4, 5: Right ascension 
and declination; columns 6, 7: K and J— K (from lKaner et al.l2006lk column 8: Ay applied to deredden the sources (see 
text); columns 9, 10: Resulting absol ute if -band magnitude and intrinsic (J—K)q; column 11: A dopted distance with i n 
betw een brackets the reference : B96:|Bronfman §\ al.| l|l996l) wit h the galactic rotation m odel oflBj^mc^^Fjli^d fl993). 



C87 
K89 



Caswell fc Hayned dl987l) C 91:IChini fc Neckell (11981^ C95:ICarpenter et al]|)l995|) . F02: lFigueredo et al.l|20Q2ji . 



vnes qiMi , T , . . . . _ _ 

Koempe et all lll989tl . K 94:lKurtz et all il9M\. K 01: lKarnik et alJ ll200ll) , Kpc: Kurtz, private communi cation, 
L92:lLiseau et al.Ml992h . S90: ISimpson fc Rubinl Jl99(il) . T91: lTapial lll99lh . WSg ^Wood fc Churchwelll lll989bh . W97: 
IWalah et al .1^1 997^. 



3.1. Position in the colour-magnitude diagram 

In Fig. the objects are plotted in a (K, J — K) colour- 
magnitude diagram. The observed if-band magnitude is 
transformed to absolute magnitude by taking into account 
the distance and interstellar reddening. The adopted dis- 
tance is listed in Tabled The solid li ne is the location of 
the zero- a ge main sequence (ZAMS) (lHanson et al.lll997l 
Ifwll2000h . 

It turns out that practically all candidate mas- 
sive YSOs are m embers of an embedded cluster (cf. 
iKaper et al1l2006f The massive YSOs located in IRAS 
06412-0105 and IRAS 15411-5352 look isolated, but those 
regions are associated on a larger scale with (massive) 
starforming regions. Our observations also do not go deep 
enough to exclude the presence of low-mass members of 
the cluster. In some cases we have been able to spec- 
troscopically ide ntify other early-type cluster members 
ijBik et alJ l2005[) . As the intrinsic colour (J - K) of 
all early-type main-sequence stars is ~ —0.2, it is then 
straightforward to derive the amount of interstellar red- 
dening. 

Assuming that the extinction does not vary a lot over 
the Hn region, this method to estimate the reddening 
could be applied to nine objects (black circles in Fig.^J. 
For the other objects (grey circles) the average (J — K) 



of the embedded cluster members was taken to measure 
the reddening, assuming that these stars are all on the 
main sequence and thus have (J — K)q ~ 0. It may well 
be that the massive YSOs are more embedded compared 
to the already emerged OB stars, and then the derived 
amount of interstellar extinction is a lower limit. However, 
the assumption that the extinction is uniform over the 
Hn region is not likely to be valid. To illustrate this with 
an example, the visual extinction measured in the direc- 
tion of NGC2024-IRS2 is about 3 m agnitudes higher than 
to NGC2024-IRS2b, just 5" away (|Lenorzer et all 120041: 
iBik et aljEofii. 

The objects in Fig. ^show a wide range in (J — K), 
from to almost 6. Also a number of objects have a very 
bright K magnitude. If we assume that the if-band flux 
can be fully attributed to photospheric emission, the abso- 
lute K -band magnitude of about half of the stars would be 
too bright to be main sequence O stars and would have 
to be classified as supergiants. The star-forming regions 
containing these objects are supposed to be young and 
unlikely to host evolved supergiants. This is an indica- 
tion that the bright if-band magnitude is due to a near- 
infrared excess. 

For most UCHII regions distance determinations are 
not very reliable. The distances are mostly based on 
the radial velocities measured in molecular lines (e.g. 
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X) 




J-K 



Fig. 1. Distribution of massive YSOs in the (K , J — K) colour magnitude diagram. The zero-age main sequence 
(ZAMS) is indicated, as well as the location of some supergiants (o). For each object an estimate of the amount of 
interstellar extinction has been obtain ed using other m embers of the embedded cluster. In some cases a spectrum of 
a spectroscopically identified OB star ijBik et al.l l2005) could be used (black circle); alternatively, an average value for 
the reddening is derived assuming that the other cluster members are main sequence stars (grey circle). The diagonal 
line in the top-right of the figure is the dereddcning line representing 10 magnitudes of visual extinction. The adopted 
distance and the uncorrected values of K and J — K are listed in Table ^ The ob jects marked by a "□" and a 
"x" symbol are th e extinction corrected positions of the massive YSOs found in M17 ' {Hanson et alJll997l) and W 31 
{Blum et al.ll200l[) . respectively. The position of NGC2024-IRS2 l|Lenorzer et alJl2004l) and S106 jFelli et alJll984h is 
indic ated by a For comparison, the ar eas occupied by Herbig AeBe stars {Eiroa et al.ll2002^ . B[e] stars {Voorsl 
1999), and Be stars l|Doughertv et al.lll994j) are indicated. 



iBronfman et al.ll996h and converte d to a distance using a 
rotation model for the galaxy (e.g. iBrand fc Blitzl ll993). 
These distance determinations give a near and a far solu- 
tion in some directions and, furthe rmore, the mod e l does 
not account for peculiar motion. IBrand fc Blitd (^993) 
show that in some directions the observed radial veloc- 
ity spread can be up to 30-40 km s _1 , which makes the 
distance determination uncertain. If a star of known spec- 
tral type is detected in the embedded cluster, another 



the 



constraint on the distance is obtained, so that, e.g. 
distance ambiguity can be resolved (cf. iBik et al1 l2i 



For 5 regions a spectroscopic parallax based on op- 
tical data has been determined. IRAS 17136-3617 and 
IRAS 18006-2422 are visible in the optical, while IRAS 
06058+2138, IRAS 06061+2151 and IRAS 08576-4334 are 
associated with an optically visible H n region or OB as- 
sociation (for references, see column 12 of table GJ. 
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Fig. 2. JT-band spectra of candidate massive YSOs in two wavelength settings: the Br7 setting (2.08 - 2.19 /Km, 
left) and the CO setting (2.29 - 2.40 fim, right). The spectra are normalised to the continuum. Note the gap in 
wavelength between the two settings. At the top of the figure the identification of the spectral lines is given. The 
CO setting displays the most prominent differences between the targets. Note the double-peaked Pfund emission lines 
in 16571nrl281. The narrow absorption lines present in the CO setting of 08576nr408 and 17258nr593 are remnant 
telluric lines. 
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Fig. 2. Normalised if -band spectra of candidate massive YSOs (continued). CO first-overtone emission is detected in 
5 objects with pronounced difference in line shape. 16164nr3636 and 18006nr766 show both CO and Pfund emission 
lines. The narrow absorption lines present in the CO setting of 07299nr598 are remnant telluric lines. 
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A (fJ,m) Element Transition 



0893 


Fe ii 


3d6.4p-3d6.4s 


1127 


He I 


ls.3p-ls.4s 


1138 


He i 


ls.3p-ls.4s 


1218 


H 2 


1-0 S(l) 


1374 


Mgn 


5s-5p 


1438 


Mgn 


5s-5p 


1624 


He I 


ls.4d-ls.7s 


1649 


He I 


ls.4f-ls.7g 


1661 


Hi(Bi'7) 


7-4 


2935 


CO 


2-0 


3227 


CO 


3-1 


3535 


CO 


4-2 


3829 


CO 


5-3 


3329 


Hi 


33-5 


3364 


Hi 


32-5 


3402 


Hi 


31-5 


3445 


Hi 


30-5 


3492 


Hi 


29-5 


3544 


Hi 


28-5 


3603 


Hi 


27-5 


3669 


Hi 


26-5 


3743 


Hi 


25-5 


3828 


Hi 


24-5 


3924 


Hi 


23-5 



Table 2. Identification of the lines detected in the K- 
band spectra of candidate massive YSOs. Column 1 lists 
the rest wavelength, column 2 the corresponding element 
and column 3 the line transition. 



3.2. Bry setting 

The objects in our sample were selected based on the pres- 
ence of a broad Br7 emission line. The spectra are dis- 
played in Fig. [21 Table [21 lists the identification of the lines 
detected in the spectra. The full width at half maximum 
(FWHM) ranges from 100 to 230 km s" 1 (the spectral 
resolution is ~ 30 km s _1 ). The large FWHM suggests 
that the emission has a circumstellar (e.g. a rotating disk) 
and not a nebular origin. In a H n region of 10,000 K the 
isothermal sound speed is 17 km s -1 and would remain un- 
resolved. We note, however, that hydrogen recombination 
lines detected at radio wavelengths in som e hyper-compact 
Hn r egions display comparable widths {Kurtz fe Franco! 
|2002|) . In two objects (16164nr3636 and 06412nrl21) He I 
emission lines are detected. In the case of 06412mT21 
these lines are not resolved, suggesting a nebular origin. In 
16164nr3636 the lines are broad and resolved like the Br7 
line, and probably produced in the circumstellar environ- 
ment. He I emission is only produced if the nebula is ex- 
pose d to the strong UV f i eld of a hot, early-type star (cf . , 
e.g., lHanson et alJ l2002t iMartm-Hernandez et alJ 12003'). 
In the majority of objects nebular, spatially extended, H 2 
1-0 S(l) emission is observed. This emission is produced 
by molecular hydrogen that is excited by either shocks or 
the UV radiation field. 

The strength of the Br7 emission covers a wide range 
(measured equivalent width (EW) from 1 to over 100 A, 



with the majority between 1 and 10 A, Table We do 
not find a relation between the Br7 strength and width, 
nor does the Br7 strength correlate with the appearance 
of other lines detected in this wavelength setting. One ob- 
ject, 06058nr221, shows a double-peaked Br7 profile, with 
a peak separation of 95 km s _1 (Fig. Such a double- 
peaked profile is also ob served (though in P a<5) in massive 
YSOs detected in M17 {Hanson et alJll997|) . The double- 
peaked nature suggests the presence of a rotating disk. 
The Br7 profile in several other targets displays an asym- 
metry. Figure |3 highlights some of the observed Br7 lines. 
Some objects show an additional red component, while 
other objects exhibit blue-shifted emission. There are no 
signatures indicating the presence of C IV and N in emis- 
sion lines, demonstrating that these stars are not "normal" 
OB super giants with a strong stellar wind causing the Br7 
emission {Hanson et al.lll996|) . 

Several objects exhibit an Fen emission line 
(2.089 /im), while a few show broad Mgn emission (2.138, 
2.144 /im). The Mgn emission is lik ely produced by the ex- 
citation through Ly/3 fluorescence llBowenlll947l) . Like the 
Fen emission, the Mgn emission lines indicate the pres- 
ence of a den se and warm (several 1000 K) circumstellar 
environment {McGregor et al.lll988|) . 

3.3. CO setting 

3.3.1. CO first-overtone emission 

One of the most remarkable spectral features found in 
these objects is the CO first-overtone emission in the 
range from 2.29 — 2.40 /Ltm. CO bandheads are detected 
in 5 sources (Fig. 0), with different shapes of the band- 
head. The objects 11097nrl218 and 08576nr292 show a 
very steep bandhead, while the other 3 objects have a 
more shallow bandhead shape and exhibit an extended 
blue wing. The shape of the bandhead reflects the veloc- 
ity dispersion of the CO molecules. In 11097nrl218 and 
08576nr292 the CO gas has a small velocity dispersion, 
which also results in a high contrast in the rotational J 
lines. In the other objects the CO molecules have a much 
larger velocity dispersion, resulting in a reduced contrast 
of the J lines and a broadening of the bandheads. 

The observed extended blue wing of the bandhead can- 
not be explained by a gaussian velocity distribution. A ke- 
plcrian velocity dis tribution, howeyer, can repro duce such 
a bandhead shape. IChandler et"aTl ill 993L liflfll have ob- 
served and modelled high resolution CO 2-0 observations 
of a sample of YSOs and modelled the kinematics of the 
CO gas. These models are based on a circumstellar accre- 
tion disk or a neutral stellar- or disk wind. Most of the 
object s display a CO bandhe ad with a blue shoulder, for 
which IChandler et alJ {l995|) obtained the best fit with 
a circumstellar disk. A detailed analysis and modelling 
of the CO first-overtone bands in our massive YSO sam- 
ple su pport this conclusion and is presented i nlBik fe Thi 
(2004!,). Similar results have been obtained bv lBlum et al. 
(2004J). Contrary to the large difference in shape of the 
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07299nr618 


4.0 


b 


b 


1.0 
















08576nr292 


67.1 


76.5 


64.5 


32.3 
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25.0 


28.7 


28.7 


19.1 
















16164nr3636 


12.4 


b 


b 


b 
















18006nr766 


17.2 


16.2 


12.2 


b 

















Table 3. Equivalent width (EW) measurements of the emission line objects. Upper table: Columns 2-7 give the EW for the spectral lines in the Br7 setting. 
The EWs are measured in Angstrom and emission is given as positive. In column 8 the FWHM of Br7 is given in km s _1 . Columns 9 and 10 show whether the 
CO and Pfund-lines are present in the spectra or not. ("-" means not present, "n" means not observed). Middle table: The EW measurements of the Pf-lines. 
The "b" means that the line could not be measured because of blending with the CO lines. The Pf-lines of 16571nrl281 are double peaked with an average 
peak separation of 189 ± 33 km s _1 . Lower table: The EW measurements of the CO first-overtone bandheads. The "b" means that the lines are blended with 
cold CO absorption lines or Pfund lines. 




Fig. 3. Differences in Br7 line morphology. The top spectrum in the left panel (06084ml 18) shows a blue extended 
wing, while the bottom spectrum in the middle panel exhibits a red extended wing of the Bi'7 profile. The bottom 
spectrum of the left panel (06058nr221) displays a double peaked Br7 line with a peak separation of 95 km s -1 , 
suggesting that it originates in a rotating disk. In the right panel, the Br7 line and Pf25 line of 16571nrl281 are 
plotted. By measuring the width of the two lines (Table |3J, Pf25 turns out to be much broader (and double peaked) 
than the Br7 line. Two telluric lines at —360 and 345 km s _1 with respect to Br7 hamper an accurate determination 
of the Br7 line wings. 



CO bandheads, the FWHM of Br7 is rather constant, 
also the CO lines with a small velocity dispersion show a 
broad FWHM of Br7 (Table EJ). This suggests that Br7 is 
not formed in the same region/geometry as the CO band- 
heads. 

In the spectra of 16164nr3636 and 07299nr618, also CO 
absorption lines are detected (2.33 — 2.35 /zm). These lines 
are the lower J-lines of the 2-0 bandhead and indicate cold 
(~ 40 K), foreground molecular gas. 

3.3.2. Pfund lines 

In five objects hydrogen Pfund lines are detected, includ- 
ing two which also have the CO first-overtone bands in 
emission. In the spectrum of 16571nrl281 the Pf-lines are 
double peaked with a separation of 189 ± 33 km s _1 sug- 
gestive of a disk origin. In Table H the FWHM of Pf25 is 
given in comparison with the Br7 line. It turns out that 
the Pf25 line is always much broader than Br7, which 
would be consistent with a disk origin (Sect. 5). 

The Pf-line emission indicates a line forming region 
of high density {N c - 10 8 cm" 3 ). The flux ratio of Pf25 
and Br7 is listed in the last column of Table 0] The ra- 
tio varies from 0.08 to 0.3, in all cases larger than what 
is expected in the optica lly thin case B scenario (< 0.02, 
IStorev fc Hummeilll995h . This means that when the emis- 
sion of the Pf25 and Br7 are emitted by the same gas, this 
gas must be (partially) optically thick. However, when the 
two lines are produced at different locations, this ratio 
provides us with the ratio of the two line forming surface 
areas assuming that the emitting material is fully optically 
thick. 

If the Pf-lines originate in a keplerian disk, they are 
likely formed in the inner disk region, while for example 



Object FWHM Br 7 FWHM Pf25 Pf25/Br 7 
(km s' 1 ) (km s' 1 ) 



08576nr292 


220 


± 


2 










11097nrl218 


191 


± 


3 










07299nr618 


138 


± 


2 










16164nr3636 


132 


± 


1 


273 


± 


8 


0.3 


18006nr766 


99 


± 


1 


282 


± 


13 


0.1 


06412nrl21 


147 


± 


1 


269 


± 


7 


0.2 


15411nrl955 


106 


± 


1 


140 


± 


37 


0.08 


16571nrl281 


204 


± 


1 


328 


± 


8 


0.2 



Table 4. The FWHM of the Br 7 line compared with that 
of the hydrogen Pf25 line. The Pf23 line is stronger than 
the Pf25, but this line is heavily contaminated by residu- 
als of atmospheric lines. For the fitting of 18006nr766 the 
Pf23 line is used, because Pf25 is not detected. In the last 
column the flux-ratio of Br 7 and Pf25 is listed. 



Br7 is formed in a much more extended part of the disk. 
In such a scenario the FWHMs of the Pf-lines and Br 7 
are expected to be different: Br 7 is predominantly formed 
in the more slowly rotating outer parts of the disk (larger 
surface area) and thus will have a smaller FWHM. 

4. Comparison to well studied objects 

In order to better understand the nature of the objects 
presented in this paper as well as the physical conditions 
under which the spectral lines are formed, the observed 
photometric and spectral properties of the candidate mas- 
sive YSOs are compared to those of well-studied objects 
in the same area of the CMD. 

In the following section we will first discuss some ob- 
jects of which the massive YSO nature has been estab- 
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lished based on a broad range of observations. After that 
we discuss the supposedly "more evolved" objects like Be 
and B[e] stars which show similar if-band spectra. In the 
last part of this section the observed properties of the can- 
didate massive YSOs are discussed in detail and compared 
to those of the reference objects. The existing spectra of 
those objects are usually of worse quality and resolution 
than the spectra of the candidate YSOs presented in this 
paper. This makes a quantitative comparison difficult. 

When objects have similar near-infrared spectra, this 
does not necessarily mean that they are in the same phase 
of stellar evolution; it does, however, indicate that both 
classes of objects provide the physical conditions under 
which a certain near-infrared spectrum is produced. 

In Table [S] the near-infrared photometric and spectro- 
scopic properties of the reference objects discussed in this 
section are summarized. Below, a short description of the 
special characteristics of the object or classes of objects is 
given. 

4.1. (Massive) Young Stellar Objects 

S106-IRS4 

S106 is one of the prototype massive YSOs. The broad 
shape of the CO bandheads in the K-band spectrum 
(Table EJ) is explained by emission from a n arrow region 
in a circumstellar disk ( Chandler et al.1119951) . The source 



is further characterised by a bi-polar outflow (v CKp = 200 
km s _1 ), with a very pronounced equatorial gap. The gap 
is also visible in the radio , where extinction does not play 
a role 1 Bally et"aTlll983|) . This indicates that the gap is 
not caused by additional extinction, but is probably the 
result of a shadow produced by a circumstellar disk, which 
is supporting the results of the fits of the CO bandheads. 

BN object 

Among the YSOs, the Becklin-Neugebauer (BN) object 
is one of the best studied. The CO bandheads are rather 
steep. If the CO emission would come from a disk, the 
modest velocity broadening leads t o the conclusion that 
the disk is viewed almost pole-on. IScoville et all <jl983h 
suggest that CO emission comes from a shock-heated thin 
layer far from the star caused by the motion of the BN 
object through a cold medium. 

NGC2024-IRS2 

The infrar e d sou rce NGC2024-IRS2 was discovered by 
iGrasdalenl Jl974h and long thought to be the ionis- 
ing sourc e of th e Flame Nebula (NGC2024). However, 
iBik et alJ l|2003l) identified NGC2024-IRS2b as the ion- 
ising source, a late O, early B star. Using simple gas and 
dust models of prescribed geometrv lLenorzer et alJ l)2004|) 
argue that the infrared flux from the circumstellar mate- 
rial of NGC2024-IRS2 is produced by a dust-free, dense 
gaseous disk. The density in the disk of NGC2024-IRS2 is 
on the order of 10 14 — 10 15 cm~ 3 in the equatorial plane at 
the stellar surface. This is 2-3 orders of magn itude higher 
than in classical Be stars l|Waters et al.lll99lj) . 

Massive YSOs in (giant) Hn regions 



A near-infrared stud y of the young Hn region M17 
by IHanson et alJ l|l997l) has revealed a sample of mas- 
sive stars possessing a near-infrared excess. The -ftf-band 
spectral properties are identical to the objects presented 
here. Spectra obtained around 1 micron revealed double 
peaked Pa<5 emission, superposed on photospheric absorp- 
tion lines. In the optical where the circumstellar material 
does not produce excess emission, the photospheric spec- 
trum of the underlying star is visible. The spectral types 
derived from these spectra are early and mid-B. 

Herbig Ae/Be stars 

The circumstellar emission of Herbig Ae stars is pro- 
duced by a passiv e , dus t y circumstellar disk (e.g. 
Mannings fc . SargenlJ Il997t iMillan-Cabet et all l200li 
Eisner et al J 120031) . For the Herbig Be stars the geome- 
try of the circumstellar material is not obvious. At sub- 
mm wavelengths the predicted emission produced by a 
cold, massive circumstellar disk is not detected, suggest- 
ing that the geometry of the circumstellar material of the 
Herbi g Be stars is differ ent from that of the Herbig Ae 
stars. lNatta et all l|2000l) suggest that this may reflect a 
difference in time scales. The dispersal of the in-falling en- 
velope by radiation pressure may trigger a rapid evolution 
of the circumstellar disk, and the disk will disappear much 
faster than in the case of Herbig Ae stars. When the Herbig 
Be stars become o ptically visible, thei r disks may already 
have disappeared. iFuente et al.l l|2003j) found evidence for 
a disk around 3 Herbig Be stars, but the derived masses 
(10 -3 — 10 -2 Mq) are substantially lower than the mass 
of the disks around Herbig Ae stars (10 -2 — 10 _1 Mq). 

4.2. Main-sequence and more evolved stars 

Be stars 

Be stars are rapidly rotating main sequence or (sub)giant 
stars that are characterised by Ha emission produced by 
high density circumstellar gas. Often the Ha line is double 
peaked. T heir infrared spec t ra are dominated by hydro - 
gen lines l|Honv et al ] l2000t IVandenbussche et alJ 12002^ . 
including many from higher energy levels. Interferometric 
observations have confirmed the hypothesis that the cir- 
cumstellar gas is distributed in a disk-like geometry 
i Vakili et alJll998UStee et alJll998i) . The physical mecha- 
nism producing and maintaining these circumstellar disks 
is not known, but is likely related to the rapid rotation of 
the central star. 

B[e] stars 

B[e] stars show hydrogen lines in emission but also include 
forbidden lines in their optical spectra, mainly of [Fe n] 
and other low ionisation species. Contrary to Be stars, 
these stars exhibit a strong near-infrared and mid-infrared 
excess. The SED peaks between 5 and 10 /im, showi ng the 
signa ture of a hot (500 - 1000 K) dust component ijVoorsI 
Il999|) . The B[e] stars comprise stars of ver y different na- 
ture. T he classification scheme proposed bv lLamers et alJ 
(199$) includes a class of supergiants as well as a class of 
pre-main sequence objects. However, for many B[e] stars 
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CO 
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K 
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0.6 
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BN 
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4. 






4- 


4- 
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NGC2024-IRS2 


BOV 


+ 


+ 




4- 


4- 


0.9 


-7.0 


G87, C93, C95, L04 
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F - late B 


+ 










1.5 


2 - -2 


E02, vdA04 


Herbig Be 


B 


+ 


? 


+ 


some 




2-6 


-1 --5 


E02, vdA04 


Be 


B 


+ 


+ 


some 




some 


0.0 - 0.6 


--3 


D94, COO 


B[e] 


B 


+ 


+ 


some 


+ 


+ 


2 - 6 f 


-6 - -lO 1 


M88,M96,V99 



Table 5. Summary of the spectral and photometric properties of the well studied objects as discussed in Sect. 21 Column 
1: Object name; column 2: Spectral type; columns 3-6: Spectral features detected in their K-band spectra; columns 
7 (J-K) value, d ereddened for foreground extinc tion; colum n 8: A bsolute, dereddened i f -ban d magnitude; column 9 : 
Refere nces: 087:lGeballe fc Persson[llil)87h C89:ICanl lll989l). C 93: IChandler et all Jl993l). C 95: IChandler et all (Il995l) 
vdA00 : lvan den Ancker et all ll2000|L S83:IScoville et alJ (|l98^) L04:lLenorzer et all ll20ol . vdA04: E02: lEiroa et alJ 
l 2002|Uvan der AnckeH<l2o7)4ji . C00:lciark fc Steele! 1)200(1) . D94: lDougbertv et al.ll|l994h . MSS ^Cregor et all (jl 98Sft . 
M96: lMorris et alJ l|l996l) . V99: IVoorsl l|l999h . t: These values are for the B[e] supergaints. 



the evolutionary status is unknown. The infrared excess of 
B [c] supergiants is likely caused by a highly non-spherical 
wind which is much dens er and slower at the equator 
than in the polar regions l)Zickgraf et alJIlQSfil). This is 
likely due to the rapid rotation of the star l|Pelur)essv et all 
l200oj) 

In Fig. H the objects discussed above are plotted, to- 
gether with the ZAMS and the massive YSOs presented 
in this paper. The Be stars are located close to the main 
sequence. The Herbig Be stars are much redder and bridge 
the gap between the Herbig Ae stars and the B[e] super- 
giants. The majority of the candidate massive YSOs are 
found near the location of the Herbig Be stars, suggest- 
ing that these objects are also of B spectral type. A small 
group of objects show a brighter if -band magnitude and 
also a bluer (J — if) colour. They are located near the 
well-known massive YSO S106-IRS4. Also NGC2024-IRS2 
is found in that area of the CMD. The location of this 
group suggests that the central stars of these objects are 
of O spectral type, i.e. hotter and more luminous than 
the Herbig Be stars. For some objects, this suggestion is 
supported by additional arguments (see below). 

4.3. Classification of our objects 

A detailed description of the fields and stellar populations 
hosting the m assive YSO candidat es discussed in this pa- 
per is given in iKaoer et alJ l|2006h . Here we focus on the 
characteristics relevant for a better understanding of the 
physical nature of these objects. In Table HO a summary is 
given of the properties of the objects. The objects will be 
discussed in groups based on their position in the CMD 

(Fig.nj. 

4.3.1. Objects with B spectral type 

The majority of the candidate massive YSOs have po- 
sitions in the CMD which suggest that they are of B 



spectral type. A few objects hardly show any infrared ex- 
cess and have colours similar to Be stars (07299nr598 and 
08576nr292, Fig. The spectral types of these objects 
(Table IH1 column 7) are estimated by subtracting a typi- 
cal if -band excess for Be stars (^0.6 magnitude). The K- 
band spectrum of 07299nr598 also has the characteristics 
of a Be star. Underneath the Br7 emission line, a broad 
photospheric absorption compone nt is detected. This i s 
also seen in the late-type Be stars l|Clark fc Steele! 2000). 
The spectrum of 08576nr292, however, does not have the 
spectral characteristics of Be stars and shows strong CO 
bandheads. 

The other objects with a larger J — K are located in 
the same area as the Herbig Be stars. None of the stars 
show similarities with the Herbig Ae stars, indicating that 
the stars we have detected are more massive. The spectral 
properties, however, are less homogeneous. Some objects 
show CO bandheads, some Pf-lines, while other objects 
only show Br7 emission. This difference in spectral signa- 
ture is, however, also seen in Herbig Be stars. 

4.3.2. Objects with O spectral type 

A small group (5 objects) is located close to S106 in the 
CMD. Their luminosity suggests that the central stars are 
O stars. They have a relatively blue colour (J — if ~ — 2), 
suggesting a similar structure of the circumstellar material 
(Sect. E3>- Mi d-infrared data of 1 1097nr693 indicate a 
very blue SED ijBarbosa et al.ll2003j) , likely dominated by 
circumstellar gas instead of dust. 

The spectral properties, however, are not similar to 
S106, but also differ among each other. Apart from the He I 
emission, the spectral features do not depend strongly on 
the luminosity of the stars, but are much more determined 
by the physical conditions, like temperature and density 
in the circumstellar material. The main effect which seems 
to depend on the luminosity is the relatively blue colour 
of the more massive objects. 
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The object which resembles the spectroscopic proper- 
ties of S106 most is 16164nr3636. The presence of cold 
CO absorption lines and the position in the CMD suggest 
that this object is deeply embedded. The extinction of the 
neighbouring OB stars used to deredden 16164nr3636 is 
obviously much less than the real extinction in this line 
of sight. The CO-bandheads in this object posses a blue 
wing which is natural ly explained by a keplerian rotating 
disk ijBik fc Thill2f)f)4]) . 

4.3.3. Objects showing strong Br7 emission 

In the majority of the candidate massive YSOs in our 
sample the Bi'7 EW is around 10 A or less. In three ob- 
jects, however, the EW is significantly higher: 06058nr221, 
06061nr676, 06084nrll8 (30-100 A, Table EJ). These ob- 
jects are all located in the lower part of the CMD (FigQ]). 
They are of equal luminosity as the Herbig Ae stars, but 
significantly redder (J — K ~ 2). 

Mid-infrared observations show that 06084nrl l8 has 
a very red mid- infrared SED l|Persi fc Tapiall2003j) . sug- 
gesting that the source is still heavily embedded. The IR- 
luminosity they derive is 350 Lq, which corresponds to 
an embedded star with a spectral type later than B3. The 
radio spectral energy distribution ha s a negative slope, is 
variable and has a non-thermal origin. iGomez et al. I ll2002l) 
suggest that the radio emission is originating in an active 
magnetosphere, like observed towards T-Tauri stars. The 
derived IR-luminosity as well as the position in our CMD 
suggest that this object is much more luminous than a T- 
Tauri star. It has a if-band magnitude similar to Herbig- 
Ae stars. 

The Br7 line of 06058nr221 is double peaked with a 
separation of 95 km s _1 , suggesting an origin in a ro- 
tating disk. The disk hypothesis has not been confirmed 
by CO or Pfund lines, as no observations were performed 
in this setting. 06061nr676 is one of the few objects be- 
ing the near-infrared counterpart to a radio UCHII region 
(see Sect. I4.4|l . which might be contributing to the Br7 
emission. 

4.4. Ionising sources of UCHII regions 

iKaper et all fcOOrt) find a near-infrared counterpart of 
an UCHII radio source in 44 % of the fields containing 
such a source . This is in agreement with the results of 
lHanson et alJ l|2002() , who show that only half of the ra- 
dio UCHII regions have a near infrared counterpart. 

The UCHII regions without a near-infrared counter- 
part are apparently still so heavily embedded that they 
are not detectable at near-infrared wavelengths. By com- 
paring the positions of the candidate massive YSOs with 
those of the UCHII regions, we find that 3 massive YSOs in 
our sample are coincident with a radio source (see Table® 
column 5), although the radio counterpart of one of them, 
06084nrll8, does not have the radio properties of a typical 
UCHII region (see Sect. ECU . 



To illustrate this, in panel (a) and (b) of Fig. 0] two 
sources are displayed with a radio counterpart, while panel 
(c) shows an example of a region where the radio source 
is not the counterpart of the massive YSO, but of a more 
embedded source. In panel (a), the K-b&nd image of IRAS 
06061+2151 demonstrates that 06061nr676 is part of a 
larger cluster and the only source coincident with a radio 
source (source A: G188.796+01.030). Panel (b) shows the 
if -band image of IRAS 06412-0105 where 06412nrl21 falls 
in the center of the radio emission. Because of the large 
distance to this source (7.3 kpc), the extent of the region 
(0.8 pc) suggests that this region is not an UCHII but a 
more evolved H 11 region. 

5. Discussion 

We have observed a large sample of candidate massive 
YSOs located in high-mass star-forming regions. The ob- 
served photometric and spectral characteristics show that 
these objects are surrounded by circumstellar material. 
The next question to address is whether this material is 
located in a circumstellar disk and, if this is the case, 
whether these disks are remnants of the accretion pro- 
cess, or formed as a consequence of the mass-loss process 
and rapid rotation such as in Be and B[e] stars. 

5.1. Origin of the spectral lines 

The different spectral lines observed in the if-band spec- 
tra of the candidate massive YSOs are formed in regions 
with different physical conditions. These conditions pro- 
vide information on the origin of the material. In some of 
the objects the spectral lines indicate that they are formed 
in a rotating disk (see Table® column 6). 

5.1.1. High-density indicators 

The special conditions required to emit the CO first- 
overtone emission are suggestive of an origin in a circum- 
stellar disk. The CO is only emitted at high density be- 
cause the critical density of the CO overtone transitions 
is high (10 10 — 10 11 cm" 3 ). The excitation temperature 
has to be high as well (2000 - 4000 K). The disk hypoth- 
esis is supported by the shape of the CO bandheads of 
two objects (16164nr3636 and 18006nr766). They show a 
CO bandhead profile displaying a blue wing, which can 
be ex plained by CO gas in keplerian rotation ((Bik fc Thi 

E3). 

CO bandheads can also be produced in a shock front. 
This i s suggested for the BN object where IScoville et alJ 
(1983) propose that the CO is formed in the shock-front 
between the Hn region and the ambient molecular cloud. 
Another tracer of shocks is H2, which is also observed 
towards our candidate massive YSOs. However, the H2 
is spatially extended, and possibly formed in the photo- 
dissociation region by UV excitation. The CO is only ob- 
served towards the point sources, suggesting that the CO 
and H2 emission are not formed by the same process. 
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Fig. 4. if -ba nd images overlaye d with radio continuum data of 3 UCH II regions. Panel ( a): if-band image of IRAS 
06061+2151 jKaper et al.l l2006t) overlayed with a VLA 3.6 cm image (|Kurtz et a.l.lll994^ . A (G188. 796+01.030), B 
(G188.793+01.030) and C (G188. 794+01. 031) correspond to 3 UCHII regions, w ith UCHII region A coinciding with 
the 06061nr676. Panel (b) if -band image of IRAS 06412-0105 l)Kaper et al.ll2006j) overlayed with a VLA 1.3 cm image 
(Kurtz, priv. comm .). Panel (c) VLT if-band acquisition image of IRAS 07299-1651, where the position of the radio 
source is indicated ijWalsh et alJll998|) . None of the three sources are coinciding with the radio source, which may be 
even more embedded. 



Pfund-lines probe high-density material which is also 
an indication for a disk-like geometry. This is confirmed 
for the Pf-lines in 16571nrl281, which are double peaked, 
indicating that the lines are emitted in a disk. The Pf-lincs 
in 16164nr3636 and 18006nr766 have a width comparable 
to that of the CO bandheads. This suggests that the lines 
share the same kinematics. They cannot be formed in the 
same region, as the Pf-lines are formed in ionised gas with 
a temperature too high to coexist with CO molecules. This 
leads to the suggestion that, in case of a circumstellar disk, 
the CO is formed in the neutral mid-plane, and the Pf-lines 
originate in the ionised upper layers. 

Two of our reference objects, S106 and NGC2024- 
IRS2, show similar bandhead profiles. In these objects also 
other evidence for a circumstellar disk is found. In the case 
of S106, evidence for a disk comes from the observed ge- 
om etry of the extende d optical and near-infrared emission 
(cf. lSmith et"aTll200l|) . For NG C 2024-IRS2, the infr ared 
SED points to a (gaseous) disk i Lenorzer et alJEooih . 



5.1.2. Origin of Br7 emission 

The majority of the objects, however, do not show Pfund 
or CO emission and therefore might have another geome- 
try of the circumstellar material. The Br7 line, however, 
which is present in all the spectra, shows a remarkably 
similar behaviour in all objects. This might suggest that 
the formation mechanism of the Br7 line could be simi- 
lar for all objects. Especially the FWHM of Br7 is rather 
constant, it varies between 100 and 220 km s . The EW 
of Bi'7 ranges from 1 to 10 A, except in 3 objects where 
the EW is much larger (upto 100 A, see Table EJ- Similar 
EW values are found in Be stars, where the larger EWs 



are found for the objects with the earlier spectral type 
(B0-B2) llClark fc Steelell2000h . 

If the circumstellar matter is, like the CO and Pf-lines 
suggest, in the shape of a circumstellar disk, Br7 emission 
may be produced in the ionised upper layers of the disk. 
The Br7 emission will be formed over a fairly large area of 
the circumstellar disk and the emission will be dominated 
by the largest surface area. This in contrast to the Pf-lines 
which are only formed in the high-density inner parts of a 
disk, which makes them much broader than the Br7 line 
(Table IU). 

The width of the Br7 line suggests that the line is not 
formed in the H n region itself. The sound speed in (ultra) 
compact Hn regions is too low to account for the Br7 
line widths. For an isothermal H n region of 10 000 K, the 
sound speed is about 17 km s _1 . We note that in hyper- 
compact Hn regions broad radio recombination l ines up 
to 180 km s _1 are observed (jKurtz fc Francoi2002(l . which 
probably reflects their expansion velocity. These objects 
are heavily embedded, much more than the UCHII regions 
and will, very likely, not be visible in the near-infrared. 

The large width of the Br7 line could also be explained 
by a bipolar outflow. The red-shifted and the blue-shifted 
lobes tilted in the line of sight would then be responsi- 
ble for the observed width. However, the other spectral 
lines which then also should be produced in these out- 
flows have different FWHMs, some are broader and some 
are even unresolved. This would be unlikely in the case of 
a bipolar outflow. Additionally, our objects show a very 
strong continuum which is not observed in that case. 

Also a normal stellar wind as observed for OB stars 
would not explain the Br7 emission lines. Only in su- 
pergiants the Br7 line is in emission. But these stellar 
winds have terminal velocities of the order of 1500 km s -1 . 
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Supergiants are not expected in young star-forming re- 
gions. If the lines are formed i n a stellar wind, Br 7 would 
be broader than the Pf-lines l|Lamers et alJll99^1 . which 
is not seen in the data (Table @J. 

It is also possible tha t the Br7 line is formed in a disk- 
wind ijDrew et alJll998j) . If the central star is emitting a 
lot of UV radiation, it will ionise the surface of the disk 
so that the disk evaporates (see Sect. 15 .3|) . This will then 
result in a dense and low velocity disk wind, where the 
outflow velocity will soon domi nate the rotation al velocity. 
The line widths predicted bv lDrew et alJ l|1998(l are of the 
order of 200 km s -1 . This scenario would also explain the 
difference in width of the CO and Br7 lines in 08576nr292 
and 11097nrl218, where the CO can be formed in a disk 
seen face on and Br7 is originating from a disk-wind. 

Summarising, the broad CO emission as well as the 
double-peaked Pf-lines strongly suggest that the material 
in these objects is located in a circumstellar disk. By com- 
paring the behaviour of Br7 in all the objects, the origin 
in a disk-wind seems the most plausible explanation. This 
means that all the objects may be surrounded by a circum- 
stellar disk. The absence of CO bandheads in the spectra 
of most of the objects might therefore be related to the 
special conditions required to emit the first-overtone emis- 
sion, rather than the abundance of the CO molecules. 



5.2. Nature of the objects 

The aim of this project is to identify the ionising source(s) 
of UCHII regions, and to determine their physical na- 
ture. The regions were selected based on the detection of 
an UCHII radio source and its connection with an IRAS 
poin t source fulfilling th e UCH II colour-colour criterion 
of IWood fc Chnrchwelll lll989al) . The near-infrared sur- 
veys ofoUCjTIIregfons performed bv lHanson et al 1 ll2002l) 
and lKaper et al.1 |2006|) demonstrate that only half of the 
UCHII radio sources show a near-infrared counterpart. 
After comparison of the positions of the candidate mas- 
sive YSOs with those of the UCHII radio positions, only a 
few objects turn out to be the actual near-infrared coun- 
terpart of an ultra-compact radio source (see Sect. 14. 3[) . 
All objects, however, are part of more extended sites of 
(massive) star formation which also host UCHII regions. 

The fact that the candidate massive YSOs are located 
inside regions of star formation suggests that these ob- 
jects are young. However, also around main sequence OB 
stars, disks are observed (Be and B[e] stars). The den- 
sities derived for the disks of Be s tars are too low t o 
emit the CO bandheads (~ 10 10 cm- 3 . IWaters et al.ll99l|) . 
Also disk-winds are not observed in Be stars. Studies 
of Be star abundances in open clusters suggest that the 
Be phenomenom mainly oc curs in the second half o f 
their main sequence lifetime ijFabreeat k. Torr ei6nll20 00l). 
InBJe] stars, CO emission is observed; iKraus fc Lamera 
(2003) show that in a very dense, non-spherical stellar 
wind, such densities can be reached. Only supergiants have 
these strong stellar winds. None of the candidate massive 



YSOs is found near the location of the B [e] supergiants in 
the CMD (Fig.UJ. 

The arguments presented above suggest that the ma- 
jority of these objects are young and their circumstel- 
lar disks are likely a product of the formation process of 
these massive stars. These disks might be the remnant 
of the extended accretion disks det ected in earlier phases 
of massive star for mation ( e.g. Shephe rd et alJ l200li 

izanol Il999t IChini et alJ 



Beltran et alJ Eool f 



20041 l.liang et alJ 



04 ICarav fc Lizauol llfli 
2005; Patel et a ,l]l200M. 



5.3. Disruption of disks around young massive stars 

In some of the star-forming regions where the mas- 
sive YSOs are detected, also massive stars are found 
which already cleare d their circumstellar environment. 
Observations of M17 llHanson et al.lll997(l as well as gi- 
ant H 11 regions l|Blumll2003|) give similar results. 

In our spectroscopic sample we find that most of the O 
stars "already" cleared out their environment while a rela- 
tively larger fraction of the B stars are still surrounded by 
circumstellar material. This leads to the suggestion that 
the most massive stars clear out their environment fastest. 

This trend is also seen when inspecting the location of 
the massive YSOs in the CMD. The majority of the mas- 
sive YSOs have a position in the CMD indicating spec- 
tral type B. A few objects, however, show evidence that 
the central star is a hot O star (see i!4.3fl . The location of 
these sources in the CMD is near S1 06 for which a spectral 
type 06-08 has been proposed (e.g. Ivan den Ancker et ail 
2000). Note that the position of NGC2024-IRS2 also sug- 
gests an O star as the central source. However, when tak- 
ing into account a strong i-T-band excess and the identifica- 
tion of the main ionising source of NGC2024, NGC2024- 
IRS2b (08) it is more likely that NGC2024-IRS 2 is an 
early B star l|Bik et al.ll2003l iLenorzer et al.ll2004l) . 

The objects where the central star likely is an O star 
show a (J — K) colour which is on average much bluer 
than those of the objects which are of spectral type B. 
This difference in colour suggests that the disk properties 
of these objects are different. O stars produce much more 
EUV radiation and possess stronger stellar winds than B 
stars, which shortens the destruction timescale of disks 
around O stars. 

iHollenbach et alJ l(l994l l2000h investigated the differ- 
ent destruction mechanisms of disks around young stars. 
For OB stars, the outer regions of the disk will be photo- 
evaporated. The extreme-UV photons ionise the top-layer 
of the disk. This gas will be heated to 10 4 K and ob- 
tain a velocity dispersion given by the local sound speed 
(~ 10 km s _1 ). Near the star this velocity will be smaller 
than the escape velocity. The gravitational radius (r g ), the 
radius at which the sound-speed is equal to the escape ve- 
locity ranges from ~ 500 AU for 03 stars to ~ 100 AU 
for early B (Tab.EJ). Beyond the gravitational radius, the 
sound speed will exceed the escape velocity so that the 
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matter can freely escape from the disk, e.g. in the form of 
a disk wind. 

This model discriminates between two scenarios, the 
weak wind case and the strong wind case. For massive 
stars, the strong-wind case is applicable. In this strong- 
wind case, the disk is flattened by the stellar wind. The 
strong stellar wind prohibits that ionised gas located be- 
yond r g is freely escaping from the disk surface, but is 
dragged over the surface of the disk to the point where 
the ram pressure of the wind balances the thermal pres- 
sure of the ionised hydrogen. From there on, the gas can 
freely escape from the disk. This strong- wind case leads to 
a faster destruction of the disk than the weak-wind case. 

The timescale as derived by iHollenbach et all l)l994|) 
depends on the number of ionising continuum photons, 
the mass-loss rate, the wind velocity and t he mass of the 
disk. T he st ellar parameters are ta ken from lMokiem et all 
l|2004 and iLenorzer et ail t004\. By assuming a disk 
mass of 2 solar m asses (see IHollenbach et all Il994t 
lYorke fc Wel3 Il99fij) . the destruction timescale ranges 
from somewhat less than 10 5 year for a disk around an 
early O star, to ~ 5 x 10 5 year for a disk surrounding a 
BO V star (Table [TJ). The assumption in this model is that 
the material in the disk is not replenished, and that ac- 
cretion has stopped. These short timescales indicate that 
already very soon after the accretion stops, the outer disk 
regions will be rapidly destroyed by the EUV photons of 
the star. 

For the destruction of the regions located within the 
gravitational radius, the direct interaction with the stellar 
wind is much more important. By considering the momen- 
tum tran sferred by the col l ision o f a spherical wind with 
the disk Holl enbach et al ] ll2f)00h derive a timescale for 
this process. Typical timescales for the destruction of the 
inner parts of the disk by the interaction with the stellar 
wind are given in column 6 of Table The mass of the 
disk inside this gravitational radius is taken to be 10% of 
the total disk mass. 

Comparison of the two timescales (TableJTJ shows that 
first the outer layers of the disk are destroyed by photo- 
evaporation and that the inner parts of the disk are re- 
moved at a later stage. So observationally one should look 
for indicators of a remnant disk by probing material near 
the star, i.e. warm gas and dust, and not the cold material 
probed at mm wavelengths. 

The detected number of massive YSOs which harbour 
an O star compared to the amount of massive YSOs which 
are of spectral type B is qualitatively consistent with the 
trend suggested by this model. The destruction of the 
disks around young O stars is much faster than those 
around B stars, so the probability to detect an O star 
surrounded by a remnant accretion disk is much smaller 
than such a disk around a young B star. This suggests that 
these objects are really young massive stars surrounded by 
a remnant accretion disk. Observations of Herbig Be stars 
also show the effect of the outer parts of the disk being 
photoevaporated by the UV photons. The disk masses, if 
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Table 7. Two different effects are responsible for the de- 
struction of circumstellar disks around O stars; the photo- 
evaporation process and the direct interaction between 
the wind and the disk. In the first 3 columns, the stel- 
lar parame ters are given used in the calculation of the 
timescales l|Mokiem et al.ll2004l ILenorzer et al.l [2004). In 
column 4, the gravitational radius is given. Beyond this ra- 
dius matter is removed by photoevaporation, while inside 
of it, the stellar wind directly removes material from the 
disk. Col. 5 and 6 give the typical times cales for these two 
processes. T he formulae are taken from IHollenbach et alJ 
jll994U200fft . 



detected at all, a re substan tially lo wer than those aro und 
Herbig Ae stars llNatta et al.ll200nHFnerite et al .1120031) . 

6. Conclusions 

In this paper, high- resolution if-band spectra of candidate 
massive YSOs located in high-mass star-forming regions 
are presented. The results can be summarised as follows: 

• The near-infrared photometric properties of the mas- 
sive YSO candidates suggest the presence of a near- 
infrared excess. This implies that large amounts of cir- 
cumstellar matter are present around these objects. 

• This prediction is confirmed by VLT/ISAAC if-band 
spectroscopy. The objects were selected on the pres- 
ence of Br7 emission, but also show other emission 
lines formed in a warm circumstellar environment like 
Fen, Mgn, Pf 23-5 - Pf 33-5 and CO first-overtone 
bandhead emission. In some objects also He I emission 
is detected. 

• In general these candidate massive YSOs are not the 
ionising sources of UCHII radio sources. Only 3 sources 
are associated with an UCHII radio source. All the 
objects, however, are located in more extended Hll re- 
gions; in most cases also OB stars are found which 
already cleared out their environment. The UCHII ra- 
dio sources located in these fields are not visible in the 
near-infrared, likely because they are deeply embed- 
ded. 

• The location of these objects inside high-mass star- 
forming regions, with associated UCHII radio sources, 
suggests that these objects are still in an early evo- 
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lutionary phase, although the objects may be more 
evolved than the ionising sources of the UCHII regions 
which often remain undetected. The massive YSOs 
found in M17 and a sample of giant H n regions suggest 
that these objects have an age of about 10 6 year. 
Most of the stars have near-infrared colours similar 
to the Herbig Be stars, suggesting that also most of 
the stars in our sample are of spectral type B. For 
5 objects, however, the location in the CMD (Fig. ^ 
suggests that the central stars are of spectral type O. 
For two sources this is supported by circumstantial 
evidence provided by the nebular emission lines and 
IRAS infrared flux. 

Photo-eva poration and disk destruction models pre- 
sented by Irlollenbach et all l)l994l l2000|) predict that 
O stars clear their circumstellar environment very 
soon after the accretion has stopped. For B stars the 
timescales are substantially longer. This is consistent 
with the detection of only a few massive YSOs with 
spectral type O. 

The Br7 line likely originates in an ionised disk wind, 
while the Pf lines are probably emitted in the inner 
regions of a disk. The CO emission lines are only pro- 
duced in a very dense medium, and are likely forming 
in the dense mid plane of the disk. 
The high density to emit the CO bandheads as well 
as the origin of Br7 in a disk-wind do not point to 
an origin in a mass-loss disk such as present around 
Be stars. These conditions are met in the dense non- 
spherical winds of B[e] supergiants, but the location of 
the sources in high-mass star-forming regions as well as 
their location in the CMD (Fig. Q are not consistent 
with a possible supergiant nature. 
The evidence presented in this paper suggests that 
the infrared excess is created by a remnant accretion 
disk, rather than a mass-loss disk. These disks are 
likely the remnants of the large ac cretion disks found - 
around high mass stars in hot cores ( Minieret^alJl^S; 
Shepherd et~aLll200lt Ir3eltran et aljl2004l IChini et all 



2004). 



Acknowledgements. AB acknowledges financial support from 
the DFG during a two-month visit at ESO Headquaters. LK 
acknowledges the support of a fellowship of the Royal Academy 
of Arts and Sciences in the Netherlands. The authors thank 
the VLT staff for support and help with the observations, 
Margaret Hanson for help with the observations and datare- 
duction process, Annique Lenorzer, Wing-Fai Thi, Stan Kurtz 
and Fernando Comeron for helpful discussions. Andrew Walsh 
and Ed Churchwell are thanked for providing the radio im- 
ages. We thank the referee for his/her suggestions to improve 
the paper. NSO/Kitt Peak FTS data used here were produced 
by NSF/NOAO. 



References 

Bally, J., Sncll, R. L. 
154 



Predmore, R. 1983, ApJ, 272, 



Barbosa, C. L., Damineli, A., Blum, R. D., & Conti, P. S. 

2003, AJ, 126, 2411 
Beltran, M. T., Cesaroni, R., Neri, R., et al. 2004, ApJ, 

601, L187 

Bik, A., Kaper, L., Hanson, M. M., & Smits, M. 2005, 

A&A, 440, 121 
Bik, A., Lenorzer, A., Kaper, L., et al. 2003, A&A, 404, 

249 

Bik, A. & Thi, W. F. 2004, A&A, 427, L13 

Blum, R. D. 2003, in IAU Symposium 212, A massive star 

odyssey: from main sequence to supernova., 458 
Blum, R. D., Barbosa, C. L., Damineli, A., Conti, P. S., 

& Ridgway, S. 2004, ApJ, 617, 1167 
Blum, R. D., Conti, P. S., & Damineli, A. 2000, AJ, 119, 

1860 

Blum, R. D., Damineli, A., & Conti, P. S. 1999, AJ, 117, 
1392 

— . 2001, AJ, 121, 3149 

Bonnell, I. A., Bate, M. R., & Zinnecker, H. 1998, 

MNRAS, 298, 93 
Bowen, I. S. 1947, PASP, 59, 196 
Brand, J. & Blitz, L. 1993, A&A, 275, 67 
Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 

81 

Carpenter, J. M., Snell, R. L., & Schloerb, F. P. 1995, 

ApJ, 450, 201 
Carr, J. S. 1989, ApJ, 345, 522 

Caswell, J. L. & Haynes, R. F. 1987, A&A, 171, 261 
Chandler, C. J., Carlstrom, J. E., & Scoville, N. Z. 1995, 
ApJ, 446, 793 

Chandler, C. J., Carlstrom, J. E., Scoville, N. Z., Dent, 

W. R. F., & Geballe, T. R. 1993, ApJ, 412, L71 
Chini, R., Hoffmeister, V., Kimeswenger, S., et al. 2004, 

Nature, 429, 155 
Chini, R. & Neckel, T. 1981, A&A, 102, 171 
Churchwell, E. 1991, in NATO ASIC Proc. 342: The 

Physics of Star Formation and Early Stellar Evolution, 

221 

Churchwell, E. 2002, ARA&A, 40, 27 

Clark, J. S. & Steele, I. A. 2000, A&AS, 141, 65 

Conti, P. S. & Blum, R. D. 2002, ApJ, 564, 827 

Cox, A. N. 2000, Allen's astrophysical quantities (Allen's 

astrophysical quantities, 4th ed. Publisher: New 

York: AIP Press; Springer, 2000. Edited by Arthur 

N. Cox. ISBN: 0387987460) 
Dougherty, S. M., Waters, L. B. F. M., Burki, C, et al. 

1994, A&A, 290, 609 
Drew, J. E., Proga, D., & Stone, J. M. 1998, MNRAS, 

296, L6 

Eiroa, O, Oudmaijer, R. D., Davies, J. K., et al. 2002, 

A&A, 384, 1038 
Eisner, J. A., Lane, B. F., Akeson, R. L., Hillenbrand, 

L. A., & Sargent, A. I. 2003, ApJ, 588, 360 
Fabregat, J. & Torrejon, J. M. 2000, A&A, 357, 451 
Felli, M., Massi, M., Staude, H. J., et al. 1984, A&A, 135, 

261 

Figucredo, E., Blum, R. D., Damineli, A., & Conti, P. S. 
2002, AJ, 124, 2739 



18 



A. Bik et al.: Massive YSOs in (ultra-)compact Hn regions 



— . 2005, AJ, 129, 1523 

Fuentc, A., Rodriguez-Franco, A., Testi, L., et al. 2003, 
ApJ, 598, L39 

Gomez, Y., Rodriguez, L. F., & Garay, G. 2002, ApJ, 571, 
901 

Garay, G. & Lizano, S. 1999, PASP, 111, 1049 
Geballe, T. R. & Persson, S. E. 1987, ApJ, 312, 297 
Grasdalen, G. L. 1974, ApJ, 193, 373 
Hanson, M. M. 1998, in ASP Conf. Ser. 131: Properties of 

Hot Luminous Stars, 1 
Hanson, M. M., Conti, P. S., & Rieke, M. J. 1996, ApJS, 

107, 281 

Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 
489, 698 

Hanson, M. M., Luhman, K. L., & Rieke, G. H. 2002, 
ApJS, 138, 35 

Hollenbach, D., Johnstone, D., Lizano, S., & Shu, F. 1994, 
ApJ, 428, 654 

Hollenbach, D. J., Yorkc, H. W., & Johnstone, D. 2000, 

Protostars and Planets IV, 401 
Hony, S., Waters, L. B. F. M., Zaal, P. A., et al. 2000, 

A&A, 355, 187 
Jiang, Z., Tamura, M., Fukagawa, M., et al. 2005, Nature, 

437, 112 

Kaper, L., Bik, A., Comcron, F., & Hanson, M. M. 2006, 

to be submitted to A&A. 
Karnik, A. D., Ghosh, S. K., Rengarajan, T. N., & Verma, 

R. P. 2001, MNRAS, 326, 293 
Kocmpe, C., Baudry, A., Joncas, G., & Wouterloot, 

J. G. A. 1989, A&A, 221, 295 
Kraus, M. & Lamers, H. J. G. L. M. 2003, A&A, 405, 165 
Kurtz, S., Churchwell, E., & Wood, D. O. S. 1994, ApJS, 

91, 659 

Kurtz, S. & Franco, J. 2002, in Revista Mcxicana de 

Astronomia y Astrofisica Conference Series, 16 
Lamers, H. J. G. L. M., Najarro, F., Kudritzki, R. P., et al. 

1996, A&A, 315, L229 
Lamers, H. J. G. L. M., Zickgraf, F., de Winter, D., 

Houziaux, L., & Zorec, J. 1998, A&A, 340, 117 
Lenorzer, A., Bik, A., de Koter, A., et al. 2004, A&A, 414, 

245 

Liseau, R., Lorenzetti, D., Nisini, B., Spinoglio, L., & 

Moncti, A. 1992, A&A, 265, 577 
Mannings, V. & Sargent, A. I. 1997, ApJ, 490, 792 
Martin- Hernandez, N. L., Bik, A., Kaper, L., Tielens, 

A. G. G. M., & Hanson, M. M. 2003, A&A, 405, 175 
McGregor, P. J., Hyland, A. R., & Hillier, D. J. 1988, ApJ, 

324, 1071 

Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, 
ApJ, 546, 358 

Minier, V., Booth, R. S., & Conway, J. E. 1998, A&A, 
336, L5 

Mokiem, M. R., Martin-Hernandez, N. L., Lenorzer, A., 
de Koter, A., & Tielens, A. G. G. M. 2004, A&A, 419, 
319 

Morris, P. W., Ecncns, P. R. J., Hanson, M. M., Conti, 
P. S., & Blum, R. D. 1996, ApJ, 470, 597 



Natta, A., Grinin, V., & Mannings, V. 2000, Protostars 

and Planets IV, 559 
Neckel, T. & Staude, H. J. 1984, A&A, 131, 200 
Patel, N. A., Curicl, S., Sridharan, T. K., et al. 2005, 

Nature, 437, 109 
Pelupessy, I., Lamers, H. J. G. L. M., & Vink, J. S. 2000, 

A&A, 359, 695 
Persi, P. & Tapia, M. 2003, A&A, 406, 149 
Rousselot, P., Lidman, C, Cuby, J.-G., Moreels, G., & 

Monnet, G. 2000, A&A, 354, 1134 
Scovillc, N., Kleinmann, S. G., Hall, D. N. B., & Ridgway, 

S. T. 1983, ApJ, 275, 201 
Scovillc, N. Z., Hall, D. N. B., Ridgway, S. T., & 

Kleinmann, S. G. 1979, ApJ, 232, L121 
Shepherd, D. S., Claussen, M. J., & Kurtz, S. E. 2001, 

Science, 292, 1513 
Shu, F. H., Adams, F. C, & Lizano, S. 1987, ARA&A, 

25, 23 

Simpson, J. P. & Rubin, R. H. 1990, ApJ, 354, 165 
Smith, N., Jones, T. J., Gchrz, R. D., Klebc, D., & Crecch- 

Eakman, M. J. 2001, AJ, 121, 984 
Stee, P., Vakili, F., Bonncau, D., & Mourard, D. 1998, 

A&A, 332, 268 
Storey, P. J. & Hummer, D. G. 1995, MNRAS, 272, 41 
Tan, J. C. 2004, ApJ, 607, L47 

Tapia, M. 1991, Memorie della Societa Astronomica 

Italiana, 62, 731 
Vakili, F., Mourard, D., Stee, P., et al. 1998, A&A, 335, 

261 

van den Ancker, M. E., Tielens, A. G. G. M., & Wesselius, 

P. R. 2000, A&A, 358, 1035 
van der Ancker, M. E. 2004, in High Resolution 

Infrared Spectroscopy in Astronomy, ed. H. U. Kaufl, 

R. Siebenmorgen, & A. Moorwood, 309-314 
Vandcnbusschc, B., Beintema, D., de Graauw, T., et al. 

2002, A&A, 390, 1033 
Voors, R. H. M. 1999, Ph.D. Thesis, University of 

Amsterdam 

Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, 

G. 1998, MNRAS, 301, 640 
Walsh, A. J., Hyland, A. R., Robinson, G., & Burton, 

M. G. 1997, MNRAS, 291, 261 
Waters, L. B. F., Marlborough, J. M., van der Veen, 

W. E. C, Taylor, A. R., & Dougherty, S. M. 1991, A&A, 

244, 120 

Wolfire, M. G. & Cassinelli, J. P. 1987, ApJ, 319, 850 
Wood, D. O. S. & Churchwell, E. 1989a, ApJ, 340, 265 
Wood, D. O. S. & Churchwell, E. 1989b, ApJS, 69, 831 
Yorke, H. W. & Sonnhalter, C. 2002, ApJ, 569, 846 
Yorke, H. W. & Welz, A. 1996, A&A, 315, 555 
Zickgraf, F.-J., Wolf, B., Stahl, O., Lcithercr, C, & Klare, 
G. 1985, A&A, 143, 421 









A. 


Rik pf al ■ 

iJlri L Oil.. 


Massive YSOs in 


(ultra-)compact H n regions 


19 


(1) 


(2) 


(3) 


(4) 


(5) 


W 








(10) 


Object 


CO 


Me 11/ 


Pf 


UCHII 


Disk 


near-IR 


IRAS 


Other 


Cluster 






Fell 








Sp. type 


Sp. type 


Sp. type 


Properties 


06058nr227 


n 




n 






B? 


B1V 




Cluster member, B 1-B2V star 




















fonnH in Qfimp pIiiq+pp 

1WL111L1 Ill OdlllL, L LI OLvl a 


06058nr221 


n 




n 




+ 


late B 


B1V 




See 06058nr227 


06061 nrfi7fi 








4. 




B? 


B1V 


B2V (K94) 


r^lll^tpT TTIPITlVtPr R1-R*?V sfflT - 
\ — J L UolCl lllClllUvl -. 1 J 1 I J £i V O Li CXI 




















found in same cluster , central 




















source UCHII region. 


06084nrll8 


n 


-(- 


n 


4- 




mid-B 


B1.5V 


>B3 (P03) 


Cluster member (GGD 14), neg- 




















fifivp radio tilnnp f C^f) 1 ? i Rprl 

<Ju L 1 VC 1 CXLllLJ olUUL I \J \J 1 . 1 i/CVJ. 




















mid-IR SED (P03). 


06084nrll4 






n 






mid-B 


B1.5V 




Cluster member. 


06412nrl21 




-f 


-f 


4- 


+ 


mid-0 


09.5V 


08V (Kpc) 


Nebular Hei emission, ionising 


















07V (N84) 


snnrrp TT("!HTT' brip^Vitpst smirrp 

□ Lllll \J Villi , L/l ItllUjOL LJ W \_Ll \_'V^' ■ 


07299nr598 












B2V 


B2V 


late B (COO) 


Likely not related to IRAS 




















cimirpp Rp titflr 

uU 111 Ljvj , IJvj O L CXI ■ 


07299nr61 8 


_|_ 








+ 


mid-B 


B2V 




Slpa+tprpfl licVii" nf mnrp pitiKpH- 

KJVjCLu IC1 vl lltLllL Wl 111LJ1 L v_-lll l_lt.Ll 




















dpd ^niirpp^ 


07299nr43 


+ 


-f 








mid-B 


B2V 




Blended with 07299nr43. 


07299nr314 




-)- 








B? 


B2V 




Deeply embedded. 


08576nr292 


-|- 








+ 


mid-B 


B1V 




Located at edge of the cluster, 2 




















09V-B1V stars in center. 


08576nr408 












early-B 


B1V 




See 08576nr292. 


1in97nr1218 

11L7C 1 111 IljIU 


_|_ 








+ 


early-B 


06V 




Lnpfifpd in pIiiq^pp nnf main 

1 J \_J Cl L C \_1 111 VjIUijuCI , 11LJ \j 111CX111 




















ionising source, peak of radio 




















emission even more embedded. 


11097nr693 












late-0 


06V 




Blue mid-IR SED (B03), lo- 




















cated Hose to 11097nr1218 

L_*drLL_*Ll LlUoL L w 11U « V 1 111 1 ij 1 U ■ 


15411nr1Q55 




4- 


4- 




+ 


mid-0 


08V 




IVpbiilar LTp t pirn'scnon main inn- 

liCU V_LL<ArL lit. 1 V31-11100 1V_J41 , lllttlll 1LJ11 




















ising source H II region. No radio 




















observations available. 


16164nr3636 


4- 


4. 


_|_ 




+ 


mid-0 


06 V 


>04V (C87) 


TV/Tain inniQincr "sniirpp vprv pom- 

4V1CX111 IWlllOlllfcl OLJLllL^L^ VL1 y L^Wlll 




















pact H II region (3.5"); nebular 




















He I emission, embedded in mas- 




















civp pnvplnnp f TCD1 1 

□ 1VC L-ll VCIL/UL I livl /. 


1 6571 nrl 281 

A.\J\J 1 1111 l^Ol 






4- 




+ 


early-B 


B0.5V 


>09V (C87) 


r^ln^tpT TTIPTTlbpr TlP"vt fn ORV- 

VJIUDILI 111C111UL1 , 11CAL LLJ L/U V 




















B1V star. 


17136nr649 












late-0 


08V 




In center of small cluster 




















(GM24) next to IJOHTT region 

\ VJX 1 VI t: / , ULiA L UU \J1111 1 CtlWll 


17258nr593 




4- 








early-B 


08V 






18006nr766 


_|_ 




4. 




+ 


mid-B 


U9.5V 




T.npfifpd in np^rf fn flip inn- 

UU L.CX L L>L1 111 4V1LJ, llt.jV L IrVJ L11C 1LJ11 




















ising source Her 36 (07V). 


18507nr248 












late-0 


05.5V 




Located in cluster, also 05-06 




















V star detected 



Table 6. Summary of the properties of the massive YSOs. Columns 2-4: Spectral lines (— : not present, +: present, 
n: not observed); column 5: Is the source the counterpart of a radio UCHII?; column 6: Spectral evidence for a 
circumstellar disk. (+: evidence for a disk); column 7: Spectral type estimated based on the position in the CMD 
(Fig. [IJ; column 8: Spectral type based on the IRAS luminosity of the IRAS source. This is an upper limit on the 
spectral type of the massive YSO, as usually more stars contribute to the IRAS flux. Colu mn 9: Spectral typ e estimate 
based on other methods. Col umn 10: Propert ies of the surround ing cluster, taken f rom lKaper et ahl fcnorf) . The OB 
star spectra are disc ussed in iBik et alJ l|2005|) . References: K94: iKurtz et alJ l|l994l radio flux UCHII region); P03: 
Persi fc Tapial 1120031 m id-infrared luminosity); Kp c: Kurtz, private com munication (radio flux UCHII region); N84: 
Neckel fc Staudd lll984l optica l spectroscopy); C00: IClark &: Steele! ()2000t Br7 absorption, similarity with Be stars); 
C87: ICa,swell fc Havnesl iFoH singl e dish radio observat i ons an d ass uming that a single star is responsible for the 
ionisation): C02: IC6mez et all <|2002lk B03: lBarbosa et all l|2003h : K01: lKarnik et alJ l|200lh . 



